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\^ We holographically engineer a periodic lattice of localized fermionic impurities 

(N 

within a plasma medium by putting an array of probe D5-branes in the background 
^ produced by N D3-branes. Thermodynamic quantities are computed in the large 

o 

1^ N limit via the holographic dictionary. We then dope the lattice by replacing some 

^ of the D5-branes by anti-D5-branes. In the large N limit, we determine the critical 

c3 temperature below which the system dimerizes with bond ordering. Finally, we argue 

that for the special case of a square lattice our system is glassy at large but finite 
N, with the low temperature physics dominated by a huge collection of metastable 
dimerized configurations without long-range order, connected only through tunneling 
events. 
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I. INTRODUCTION 

Recently there has been a flurry of activity applying the anti-de Sitter (AdS)/conformal 
field theory (CFT) correspondence [IHS] to the study of condensed matter systems. Holo- 
graphic systems where the CFT exhibits strongly coupled avatars of metallic phases [1], 
superfiuid-insulator transitions [5j, superconductivity [S^iTj, and Fermi liquid theory [SHU] 
have all been proposed. Also proposed are gravity duals for scale-invariant nonrelativistic 
field theories, enjoying Galilean invariance fl2\ [13] or the smaller Lifshitz symmetry group 
[13]. For reviews, see [T3HT7] . One of the driving forces behind such vigorous activity is pos- 
sible applications to condensed matter systems with quantum critical points, such as heavy 
fermion materials and possibly cuprate superconductors. The AdS/CFT correspondence 
itself may even be thought of as a prime example of quantum critical phenomena where, at 
the conformal fixed point of a family of field theories, there appear emergent gravitons. At 
a more concrete level, one significant advantage the AdS/CFT correspondence offers over 
other methods of analyzing model field theories is the ability to compute real-time correla- 
tors, and hence to gain a handle on transport properties, in strongly coupled field theories. 
For a review of this aspect with relevant references, see [18] . 

One limitation of the studies to date is that the systems considered so far in the literature 
are spatially homogeneous. While such homogeneous systems might suffice for studies of 
critical phenomena, where the correlation length diverges and the microscopic structure 
of solids does not play any role, these holographic toy models are completely inadequate 
when it comes to questions involving the underlying lattice structure of condensed matter 
systems. In the first part of this paper, we attempt to remedy this situation by explicitly 
constructing holographic systems endowed with periodic lattice structures, in the context of 
type IIB string theory. 

To this end, we first consider a probe D5-brane in the AdSs x 5*^ background geometry, 
wrapping AdS2 x S'^. From the boundary field theory's point of view, this corresponds to 
adding localized fermionic degrees of freedom residing on a pointlike impurity coupled to 
A/" = 4 supersymmetric SU{N) gauge theory. This is a particular example of the more 
general structure of defect conformal field theory, investigated in the AdS/CFT context in, 
for example, [T9] - [2T] . Placing an array of such D5-branes in the asymptotically AdSs x 5*^ 
black brane geometry, we obtain a lattice of impurities immersed in the A/" = 4 plasma 
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medium at finite temperature |22]. Using the holographic dictionary, we can compute ther- 
modynamic quantities of this system in the large and large 't Hooft coupling limit, where 
the gravitational description is accurate. 

We then proceed further and dope the system, replacing half (say) of the D5-branes by 
anti-D5-branes. This doping introduces an interesting phenomenon of dimerization: whereas 
both D5- and anti-D5-branes go straight down into the black brane horizon at high temper- 
ature, at low temperature they pair up by connecting with each other far from the horizon. 
Again, in the large N and large 't Hooft coupling limit, we can compute the free energy 
and determine the critical temperature at which such dimerization occurs. With ease, we 
can also engineer systems with various plateaux of temperature ranges within which only a 
fraction of the unit cell of the lattice dimerizes. In the special case of a square lattice, we 
argue that our system exhibits glassy behavior at large but finite A^. 



The organization of the paper is as follows. In Sec II we build a holographic lattice and 



study its thermodynamic properties. Section III is an intermezzo where we consider a pair of 
a D5-brane and an anti-D5-brane in the black brane background. Equipped with the result 



in Sec, in, in Sec, IV, we engineer various systems which dimerize at low temperature, and 
determine the critical dimerization temperature (s). In SecjVj we briefly argue that a special 
case of our construction yields a glassy system. We conclude by suggesting numerous future 



directions in Sec VI Appendix A summarizes our conventions. 



II. HOLOGRAPHIC LATTICES 



In this section, we gradually work up to a lattice of probe D5-branes in the black brane 
background, while explaining its field theory dual and calculating associated thermodynamic 
quantities. The gravitational side of our holographic systems is governed by the type IIB 
supergravity action plus a probe D5 action for each D5-brane, S = Sub + J2i 'S'ds, (see 
Appendix A for details). For simplicity, we will not consider coincident D5-branes in this 
paper, and forces between D5-branes, which vanish at zero temperature as the configuration 
is Bogomol'nyi-Prasad-Sommerfield (BPS), but would be induced at finite temperature, will 
not matter at the order to which we work. 
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A. Black brane background 

The type IIB supergravity is well known to have the following asymptotically AdSs x 5*^ 
black brane solution: 

e 

B(2) = C'(o) = C(2) =0, — = 1, (2.1) 

9s 

g^.dx^dx'' = -f{r)de + jr^ + ^^(JZ^a;^) + L\de^ + sin^dn\) (2.2) 

■'^ ' i=i 



with /(r) = (^1 - and = AT^g^Na'^ 



(2.3) 



F(5) = dC(4) = -(1 + *) — ((it A dr A (ia;i A dx2 A dxs). (2.4) 



For later convenience, we foliated 5*^ by S'^s parametrized by their latitudes 6', and dVL\ is 
the metric on a unit 4-sphere. This solution can be obtained as the decoupling limit of 
the geometry generated by a stack of nonextremal D3-branes and, as such, it describes 
A/" = 4 supersymmetric SU{N) gauge theory at finite temperature, with the Yang- Mills 
(YM) coupling constant of the gauge theory given in terms of the string coupling constant 
as 9ym = ^T^9s [22]. 

The field theory's temperature is identified with the Hawking temperature of the ge- 
ometry above, namely, T(r+) = ^[-;^-t\/-9tt\\r=r+ = The free energy F[T] of the 
macroscopic configuration can be computed through the holographic dictionary, which reads 
in the classical limit as 

e-'^^m = Z[T] = e-^['^l (2.5) 

Here, I[T] is the properly dimensionally reduced (4 + l)-dimensional Euclidean action, eval- 
uated for the configuration with temperature T. This is by now a standard practice and the 
result is 

F™|r(,v)] ^ -^^V, (^) ^ -^^V,TV, (2.6) 

where V3 denotes the 3-dimensional volume in it direction. 



B. D5 embedding 



Let us now place a single D5-brane in the black brane background above. There is a 
one-parameter family of analytic solutions for static D5 embeddings, going straight down 
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into the black brane, characterized by latitudes 9y of 5*^ that they wrap. The position 1^ in 
the "af directions and the direction of the within 5*^ give other parameters as well. But 
for a single D5-brane, we can take these to be fixed at any value without loss of generality, 
whereas for multiple D5-branes to be considered below, their relative values would matter. 
Explicitly, by parametrizing the 4-spheres by (^i, ^2, ^3, 04) and denoting the coordinates on 
the D5-brane by ^ = (r, p, 6i, 62, 0^, ^4), the embeddings are given by^ 

{T,p,9i,92,93,(f)4^) {T,p,~^,9u,9i,92,93,(f)4) with pG(0,oo) and rG (-00,00), (2.7) 

where the worldvolume flux depends on 9^, and is given by 

{2na'J^)rp = —{2na'J^)pr = cos9y with all other components vanishing. (2.8) 

Actually, as explained in |23], possible values of the latitude 9^ G (0, vr) are quantized: 
defining a parameter v G (0, 1) by 

u=-{9y- sm9^cos9^) , (2.9) 

we have the quantization condition 

n = uN eZ. (2.10) 

In the dual field theory, introducing such a D5-brane in the bulk corresponds to coupling 
localized massless fermions to the A/" = 4 gauge theory as 



^xl^tx' + xl (Aoit, ))l + v^Mt, ))l , (2-11) 



'S'field theory — >SV=4 + J dt 

where SV=4 is the action for A/" = 4 supersymmetric SU{N) gauge theory [21], and 0/ 
are its gauge and scalar fields in the adjoint representation of SU{N), respectively, and is 
a unit 6-vector determined by the direction of the within wrapped by the D5-brane. 



The parameter n = uN, on the other hand, determines the number of fermions at the site 

correspond 

n\{N - n)! 



if [25]. More precisely, it corresponds to taking the ensemble with the density matrix 



bl<...<fe 



\bi,...,bn){bi,...,bn\ (2.12) 



^ We choose the orientation of the spacetime so that etri23e9i6i2*'304 — s-iid those of D5-branes so 

that erpe^e,e,<p, = +^-5(mduccd). 
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where, schematically,^ 

\h,...,K,) = xl...xl\0)- (2.13) 

Now we would like to compute the D5-brane contribution to the free energy. Turning 
it into the appropriate Euclidean configuration, the Euclidean action of the D5-brane is 
evaluated to be l26]-[28]^ 



.(bulk) L^sm^e^vol{S^) 1 

4 = (2vr)5,,a-3 f('^^^o.-r,), (2.14) 



where we regulated it with the regulator rcutos- We need to holographically renormalize this 
quantity. To this end, we subtract off the Euclidean action of the analogous D5 embedding 
into pure AdSs x spacetime. 



.(counter) _ ^ C7^V01(0 ) 1 

" (27r)5^,«'3 y,^y)l^cutoffj, l^.ii^j 

where we take its temperature to be T'(T) given by 

so that the induced metric on the surface at Tcutoff is the same as in the relevant black brane 
configuration [22] • This procedure is justified in the context of holographic renormaliza- 
tion [29] because both configurations are the solutions to the same equations of motion with 
the same leading Dirichlet boundary condition at r = rcutoff- Thus, the renormalized D5 
Euclidean action is given by 



.(renormalized) _ /.(bulk) .(countcr)N _ L siu ^^j^Vol(S' ) r+ 

D5 ~ UdS ~ -'D5 )— fo^\5^ T^' 



(2.17) 



Hence the leading D5 contribution to the free energy is 

r (renormalized) 



D5 = T/(— ^-'^'^ = - ^'f^'ff"^(f) r+ = -^AtTaT, (2.18) 



(27r)^(yfsa'3 Stt 
where we introduced the 't Hooft coupling A = Qyu^ ■ 



The states as written in Eq.(2.13) are not invariant under small gauge transformations, and they need to 



be supplemented by nonlocal operators for any nonzero coupling. 
^ To get correct results for free energies, it is very crucial to add the appropriate surface term for the D5 
action, analogous to the Gibbons-Hawking surface term for the Einstein-Hilbert action, so that we can 
consistently impose a Dirichlet boundary condition at r = f cutoff ■ 
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C. Lattice of D5-branes 



Finally, let us consider an array of D5-branes. For simplicity, we consider a cubic lattice 
with lattice spacing a, and let the D5-branes wrap S^s of the same direction in 5*^ and at the 
same latitude."^ On the boundary field theory side, we are just constructing a cubic lattice 
of impurities. Assembling the results of preceding subsections, the free energy of the system 
is given by 



F[T] = V, 

and hence the entropy is 

dF 



iv2 ( 



sin^^^ 1 



;T + ... 



iVMyT3 + ... 1 + 



XN 



sin^e^ 1 
Sir 



(2.19) 



(2.20) 



where dots in each parenthesis indicate suppressed terms coming from stringy corrections 
in the bulk language. When expressed in terms of boundary field theory's quantities speci- 



fying bulk boundary conditions, as above, Eqs.(2.19) and (2.20) capture the whole leading 



contributions in the large and large 't Hooft limit, up to order A^^. In particular, backre- 
action of D5-branes onto the geometry mingles only into the order A^'' contribution together 
with other bulk quantum effects. For a formal argument within the context of Einstein 
gravity and toy examples of how this works, see [30].^ 

The formulas above obtained through the bulk geometric analysis, however, have a lim- 
ited range of validity: at extremely low temperature (small horizon radius), our geometric 
analysis breaks down. This happens when the entropy contribution of the D5-branes is 
comparable to that of the background black brane geometry, namely, when T^a^ ~ The 
reason behind this breakdown is that the backreaction of D5-branes at the horizon of the 
black brane cannot be treated as small perturbation any more, since the density of D5-branes 



at the horizon 



, ^i^, increases as r+ decreases. Hence our formulas above are not applicable 



^ This guarantees the supersymmetry of the zero temperature configuration. 

^ The basic idea is that one expands the full solution to the coupled supergravity and probe action in a 
power series in and then substitutes this expansion back into the action to obtain the free energy. At 
order there are two potential contributions: the probe action evaluated on the leading solution and 
the correction to the contribution from the supergravity action due to the ^ correction in the background 
geometry. The latter however vanishes due to the equations of motion obeyed by the background. 
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outside the regime 

(2.21) 




Nevertheless, there is one "thermodynamic" quantity which we can evaluate at zero 
temperature, namely, the entropy. This comes from the fact that our lattice preserves 
supersymmetry at zero temperature, and getting the entropy is tantamount to counting the 
number of BPS configurations with given charges. Carrying out the analysis, the entropy 
at zero temperature is 

AT! \^] fV3\, {f N\ 



which is applicable for any and (^ym, and it asymptotes to (^) x N x {— z/lnz/ — (1 — 
i/)ln(l — z/)} for large A^.^ This should be contrasted with the finite temperature results which 
yields the leading D5 contribution to the entropy of order Ny/X, enhanced by V^- This is 
analogous to the strong coupling enhancement of fundamental matter in the deconfined 
plasma observed in as opposed to the zero temperature conformal medium. 

III. INTERMEZZO: D5/ANTI-D5 PAIR 

In this intermezzo, we consider a pair of a D5-brane and an anti-D5-brane immersed in the 
black brane background. In the bulk gravitational description, the anti-D5-brane differs from 
the D5-brane in its orientation, which manifests in a reversal of sign of its coupling to the 
background Ramond-Ramond fields. On the boundary field theory side, localized fermions 
introduced by the anti-D5-brane transform in the antifundamental representation of SU (N) 
while those introduced by the D5-brane transform in the fundamental representation.^ 

We will consider one D5-brane placed at 1?" = (+^, 0, 0) and one anti-D5-brane placed at 
(— 0, 0), both of them wrapping S'^s of the same angle and latitude as r — )■ oo. Intuitively 
(see Fig.l), the "ends" of both the D5- and anti-D5-branes are sucked into the black brane 
when the horizon radius r_|_ is large (or the separation Ax is large), whereas they reconnect 
outside the horizon for small horizon radius (or small separation). 

® The large entropy at zero temperature is often thought to be "unphysical." Here we interpret it just as a 

pecuharity of highly supersymmetric systems. 
^ We work in a convention in which the relation between the number of fermions introduced by the anti- 



D5-brane n = vN and its latitude 9^ is given by Eqs.(2.9l and (2.10) 



(a) 



(b) 



D5^ 



D5 



XT 



FIG. 1: Predimerization transition. (a)Disconnected configuration dominates at high temperature. 
(b)Connected configuration dominates at low temperature. 



A. Disconnected configuration 



The obvious candidate stable configuration of such a pair is just two separated configu- 



rations of the sort considered in Sec, II B with 6^ = 6p [see Fig.l(a)]. Its free energy is just 
twice that of the single D5-brane: 



Note that it is independent of the separation Ax. 



(3.1) 



B. Connected configuration 

Another candidate solution with the given boundary condition is a reconnecting solution 
[see Fig.l(b)]: a reconnecting D5-brane starts at r = oo with 1?" = (— 0, 0), dips into the 
bulk, and then comes back to r = oo now with it = (+^,0,0), effectively reversing its 
orientation as it should.^ Explicitly, we have 



e{p) = and {2TTaJ^)rp = cosO^ . =\ { ^] , (3.2) 




Incidentally, this is the reason why a pair of two D5-branes cannot reconnect. 
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FIG. 2: j^Ax as a function of k. 



and now 1^ is also a function of the coordinate p on the D5-brane. As r(p) starts from oo 
and goes toward the turning point rtum = {r^ + sfn^e 

f.r(p) 



^(p) = -(fc,0,0)x / dr' — 



1 



(3.3) 



while as r(p) goes back up from rtum and ends at oo, 

r.r(p) 



^(p) = +(fc,0,0)x / dr' — 

•J rtum ^ 



1 



sm t'l. - -74—4 



(3.4) 



Here the parameter k is implicitly determined through 



L2 



Ax 



2A; 



(1+A;2)3 



with /c = 



rlsin^6'. 



-k. 



(3.5) 



The function on the right-hand side is plotted as a function of k in Fig.2. Note that, for a 
fixed distance, such a solution does not exist for large enough r+. Also, for small enough 
r+, there are two distinct solutions with different k. By computing their free energies, we 
will determine which configuration among the two is dominant. 



We can follow the same regularization procedure as presented in Sec II B to get the free 
energy of the reconnecting solution, which yields 



connected 



^ L^sin'g^vol(g^) ^ 
(27r)5^,a'3 ^^ 

2 -, — — r I 



:i + p)^ + 



dz 



(l+fc2)I 



z^-l~k^ 



1 



(3.6) 



connected 



Ck). 



(27r)5(;,«'3 

The function -Fconnected(fc) is plotted in Fig. 3. Comparing with Fig.2, we see that when r+ 
is low enough to allow two solutions for a given Ax, the configuration with higher k, and 
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FIG. 3: -Fconnected(^) as a function of k. 

hence the one that goes less deeply into the bulk, is the dominant configuration among 
the two since it has smaller free energy. All the calculations in this section are virtually 
identical to the discussion of Wilson lines at finite temperature from fundamental strings in 
AdS spacetime [321 [33], the only difference being the overall prefactor which for us is set by 
the effective string tension of the wrapped D5-brane, whereas for the Wilson lines it was set 
by the fundamental string tension. 



C. Predimerization transition 



Here comes the punch line of this intermezzo: for a fixed separation Ax, there is a large 
phase transition as we decrease r_|_. To see this, let us start in the regime where ^^^^ 
is very large. In such a regime, as we noted earlier, there is no reconnecting solution and 
thus the disconnected configuration dominates [Fig.4(a)]. As we lower r+, there appear two 
reconnecting solutions, but at the onset they both have free energies higher than that of the 
disconnected solution and thus stay subdominant [Fig.4(b)]. As we lower r+ further and 
reach the value = fcc ~ 2 or 

^ - 0.7, (3.7) 

we see one of the reconnecting solutions, namely, the one which goes less deeply into the 
bulk, emerges as the most dominant configuration [Fig.4(c)]. Noting that free energies of 
the D5-branes have a multiplicative factor A^, we see that this marks a phase transition in 
the N = oo limit. 

This N = oo phase transition will be rendered to be a crossover for large but finite A^. 
We may avoid this by building a (3+l)-dimensional lattice whose unit cell has pairs of D5- 
and anti-D5-branes. We expect this system to have a genuine first-order phase transition in 
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(a) 



F 



Cb) 




configurations 



configurations 




configurations 



FIG. 4: Free energies as functionals of D5 configurations at various temperature ranges, (a) At 
high temperature, there is only one extremum corresponding to the disconnected solution. (b)At 
intermediate temperature, there appear two new local extrema corresponding to two connected 
solutions. (c)At low temperature, the stable connected solution dominates. 



the thermodynamic limit V3 — )■ 00, surviving a journey away from the N = 00 limit. 



IV. HOLOGRAPHIC DIMERS 



Equipped with the result in the previous section, we can now build holographic dimers, 
systems which dimerize at low temperature. Dimer systems are extensively studied in the 
condensed matter physics community, especially as a part of attempts in explaining cuprate 
superconductivity, and in trying to get a handle on the physics of Mott insulators 
And since the physics responsible for cuprate superconductivity is believed to be primarily 
(2+l)-dimensional (justifying the use of a Lawrence-Doniach model) [33], such studies have 
focused on (weakly interacting layers of) (2+l)-dimensional dimer systems. To keep the close 
analogy with such studies, we will focus on (2+l)-dimensional lattices, though it should be 
kept in mind that we are thinking of eventually stacking such (2+l)-dimensional lattices 
so that we can evade the Coleman-Mermin-Wagner-Hohenberg theorem and have a phase 
transition at finite temperature. 
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FIG. 5: A simple model that experiences the phase transition analogous to the Neel-VBS phase 
transition. 

A. Dimerization through bond ordering 

Let us consider the lattice (see Fig. 5) in which D5-branes are placed at 



a 



-^^+l]+ii^ + 9o)n,,2ny)] and a { ( ] + {{1 + go)n,,2ny) } (4.1) 



2 2J J 1 V 2 2 

and anti-D5-branes are placed at 

+((l + ^7o)n.,2n^)| and a | (^-|, -0 + ((1 + (7o)ri., 2n,) | (4.2) 

with {n^^ny) G 7? and < < 1. Working in the large and large 't Hooft coupling 
limit, then, as we cool down to ^'^"f" 0.7 or 

T.^— , (4.3) 
ago 

our lattice dimerizes through bubble nucleation of dimerized regions. This transition is 
analogous to the Neel-VBS (valence bond solid) phase transition [31], though here it is a 
first-order phase transition as a function of temperature, not a quantum critical second-order 
phase transition at zero temperature. Keeping ^ << 1, this transition happens away from 
the regime where backreaction becomes important, justifying our geometric analysis. 



B. Plateaux 

To illustrate how easy it is to engineer variants of the simple VBS dimerization described 
above, let us now consider the lattice (see Fig. 6) in which D5-branes are placed at 

■f ,+^) +((l + ^?2)n.,2n,)| and a[U^,-]\+{{l + g^)n^,2ny)\ (4.4) 
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■ 1 




FIG. 6: A model with a plateau of temperature range within which only half of each unit cell 
dimerizes. 



and anti-D5-branes are placed at 

+((l + ^72K,2n,)| and a | (^-|, -0 + ((1 + 2n,) | (4.5) 

with {nrc,ny) G and < gi < g2 < 1- Then, again working in the large and large 't 



Hooft coupling limit, as we cool down to Tic 



0.2 

a.gi ' 



only a half of each unit cell dimerizes, 



2c ~ the entire cell dimerizes. In this way, we can 



while as we cool further down to 
engineer a lattice with plateaux of temperature ranges within each of which only a fraction 
of each unit cell dimerizes. 



V. HOLOGRAPHIC GLASSES 



Finally, let us consider the case of the square lattice, which is a special limit of the model 



considered in Sec, IV A with qq = 1. For illustration, let us take = 10^°° and \/~X = 10. Be- 
low the critical temperature Tc ~ the system again dimerizes through bubble nucleation. 
After such dimerization, "spinon-pair excitation," corresponding to cutting a connected D5- 
brane into a pair of separated D5- and anti-D5-branes, is highly suppressed by a factor 
of e~^°^°° and thus the low temperature physics is dominated by dimerized configurations. 
But now there are infinitely many of them! The simple VBS-like configuration [Fig. 7(a)], 
a "resonance excitation" [Fig. 7(b)] of that configuration, another "resonance excitation" 
[Fig.7.(c)], and much more complicated configurations [like the one in Fig.7(d)] all have the 
same free energy at leading order in A^, starting to differ only at order N^. Furthermore, 
those configurations are connected only through highly suppressed tunneling events of brane 
recombinations in the bulk. Especially, note that even moving the "resonance excitation" 
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FIG. 7: Multitude of dimerized configurations for holographic glasses. (a)VBS-like configuration. 

(b) "Resonance excitation." (c)Another "resonance excitation." (d) Generic configuration. 

from site [Fig.7(b)] to site [Fig. 7(c)] involves tunneling, and thus such excitations should not 
be thought of as usual quasiparticle excitations. This marks a stark difference from standard 
dimer models and more closely resembles glassy systems, with the low temperature physics 
dominated by a humongous number of metastable dimerized configurations. We also believe 
that this system does not have any long-range order: interactions between dimers are pre- 
sumably mediated by light fields in the bulk, but, due to the peculiar geometric properties 
of AdS spacetime, such interactions are effectively rendered short ranged. Therefore, at 
large but finite N, we propose that the system transitions into a glassy phase at Tc ~ 
randomly selecting one of the plethora of dimerized configurations as we cool it down. 
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VI. FUTURE DIRECTIONS 

In this paper, we engineered holographic lattices, holographic dimers, and holographic 
glasses. In each instance, we chose to discuss a simple case to facilitate our analysis. However, 
by varying the dials in these models, we can explore a large landscape of possibilities, looking 
for qualitatively new phenomena. 

First, we can freely dial various parameters that naturally appear in the system. The 
number of D5-branes (anti-D5-branes) Nf (Nf), their directions , and the number of 
localized fermions uN {vN) can be chosen independently for each site within a unit cell. 
Also, one can think of giving a mass to the localized fermions by separating the D5-branes 
from the D3-branes in the radial direction, although the D5-brane configurations for such a 
case are analytically known only at zero temperature [36]. Turning on background magnetic 
fields [37] and/or chemical potentials for various t/(l)s is another option, and we expect to 
see the dimerization transition at zero temperature. We should also mention that localized 
fermions can be replaced by localized bosons by replacing D5-branes by appropriate D3- 
branes [23 ES]- The lattice structure can be chosen at will and does not even have to be 
periodic. 

Second, we can consider more general background geometries of the form AdSs x X5 where 
X5 is an arbitrary 5-dimensional Sasaki-Einstein manifold with volume of order L^, following 
the spirit of For such studies, flux quantization conditions analogous to the one found 
in [23] need to be worked out for general X5. We can also explore other dimensions by, for 
example, considering AdS4 x Xj and M2- or M5-branes in it, where X7 is a 7-dimensional 
Sasaki-Einstein manifold. Alternatively, we may also construct "phenomenological models," 
abstracting out the essential ingredients of defects in AdS spacetime to focus on the appli- 
cations of interest instead of the details of the string theory engineering, following the spirit 
of recent investigations of Fermi liquids [TU] . 

Third, we can try to engineer more elaborate systems with desired properties. In partic- 
ular, it would be nice if we could make the fermions hop between the lattice sites so that we 
might have a holographic dual of the Mott transition. A conceptually straightforward but 
practically hard way to accomplish this is to cutoff our geometry at some radius or to consider 

D3-branes and a lattice of D5-branes on M4 x CYg where CYq is a compact Calabi-Yau 
manifold, so that the D5-D5 strings are not completely decoupled. However it would be nice 
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if we can add such hoppings without considering compact Calabi-Yau spaces [ID] . Once such 
hopping is added, if we can completely liberate localized fermions away from lattice sites in 
some phases but not others, we have a chance of having a holographic dual representation 
of transitions in which the volume of the Fermi surface jumps [H]. 

Turning to technical challenges, it would be interesting if we can develop reliable calcu- 
lations of the free energy at order A^°. At this order, the computation is plagued by bulk 
quantum effects, but the interaction between dimers should also become evident. This could 
lead to a deeper understanding of our mysterious holographic glasses by, among other things, 
giving a small splitting between the different approximate ground states. 

Finally, one of the most interesting handles on ground states of strongly correlated sys- 
tems is the entanglement entropy. The elegant formulation of Ryu and Takayanagi allows its 
computation in the gravity dual of a strongly coupled field theory, for simple enough choices 
of the boundary geometry partitioning the system into two subsystems [12]. For bipartite 
lattice systems, another useful notion of entanglement entropy, the "valence bond entangle- 
ment entropy," which is a priori distinct from the von Neumann entropy, was introduced in 
|43j . It may be interesting to consider such generalized notions of entanglement entropy in 
our gravity duals of bipartite lattice systems. 



Acknowledgments 

We thank Stephen Shenker for very useful discussions. S.K. is also grateful to Allan 
Adams, Oliver De Wolfe, and Igor Klebanov for important remarks. This work was substan- 
tially initiated at the KITP miniprogram "Quantum Criticality and the AdS/CFT Corre- 
spondence" in June of 2009, and we would like to thank the organizers and participants for 
the very stimulating atmosphere they provided. S.K. also acknowledges the hospitality of 
the Aspen Center for Physics, and is supported by the NSF under Grants No. PIIY-0244728 
and No. PHY05-51164, and the DOE under Contract No. DE-AC03-76SF00515. A.K. is 
supported by the U.S. Department of Energy under Grant No. DE-FG02-96ER40956. S.Y. 
is supported by the Stanford Institute for Theoretical Physics and NSF Grant No. 0756174. 



18 



Appendix A: Conventions 



We set h = c = ks = I- 

Our conventions are basically those of but we will exclusively work in the follow- 
ing frame, akin to Einstein frame: {g^u)onvs = gie'"^{gfiu) string) (-S^;/)ours — (-S/ji^) string) 
(C(p))ours = fl's(C(p))string) and (e*)ours = (e*)string- This frame will nicely segregate out 
factors of Qs- Explicitly, in our frame, 

^ C(4) ^H ^ dC(2) (Ai) 



s- 



IIB 



2(27r)V4c/2 
with H = dBi2), -F(3) 



dC{2) - C{Q)H, Fi^r^) 



dC^i) - ^C(2) A + ^5(2) A dCi^2) (A2) 



and 5*05 = Sdbi + •S'cs where the Dirac-Born-Infeld (DBI) action is given by 



5i 



DBI 



I f d'^\ i - 



\ 



-det 



ab 



9s 



gab + 5(2)a6 + 27ra' J", 



ab 



(A3) 



and the Chern-Simons (CS) action is given by 

ys J D5 „_n OAR 



(2vr) 



(A4) 



' P=0,2,4,( 

with the understanding that gab, -8(2)) (^(p)) and e* are all induced from background geometry. 
Here, equations of motion that follow from varying the above action need to be supplemented 
by the self-duality condition on F(5), namely. 



(5) 



(A5) 



We will not explicitly write down appropriate Gibbons-Hawking surface terms. 
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